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. Abstract 



The effect of Lorentz invariance violation in cosmology has attracted a considerable amount of 
attention. By using a dynamical vector field assumed to point in the bulk direction, with Lorentz 
invariance holding on the brane, we extend the notation of Lorentz violation in four dimensions 
[28] to a five-dimensional brane-world. We obtain the general solution of the field equations in 
an exact parametric form for Bianchi type I space-time, with perfect fluid as a matter source. 
We show that the brane universe evolves from an isotropic/ anisotropic state to an isotropic de 
^ | Sitter inflationary phase at late time. The early time behavior of anisotropic brane universe is 

largely dependent on the Lorentz violating parameters /3i,i = 1,2,3 and the equation of state of 
' 55 ' the matter, while its late time behavior is independent of these parameters. 
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: 1 Introduction 

O , 

Brane-world scenarios with a 3-brane identified with the observable universe which is embedded in 
a higher-dimensional space-time provide an alternative to the standard cosmology in four-dimension 
[1, 2, 3]. A essential requirement on these models is that they should reproduce the important obser- 
vational cosmological data, the age of the universe, abundances of elements produced in primordial 
nucleosynthesis, etc. One of the most successful of such higher-dimensional models is that proposed 
by Randall and Sundrum whose the bulk has the geometry of an AdS space admitting Z2 symmetry 
[1]. They were successful in explaining the hierarchy problem: the enormous disparity between the 
strength of the fundamental forces. The Randall Sundrum (RS) model has had a large impact on our 
undersetting of the universe and has brought higher-dimensional gravitational theories to the fore. 
In RS type models, all matter and gauge fields live on the brane, while gravity can propagate into 
the bulk. Using the Israel junction conditions [4] and the Gauss-Codazzi equations, one can obtain 
the gravitational field equations on the brane, as employed by Shiromizu, Maeda and Sasaki (SMS) 
[5]. There are two important difference that result from the effective four-dimensional gravitational 
equations on the brane. The first one is quadratic energy-momentum tensor, TTav, which is important 
in a high-energy universe and the second one is the projected Weyl tensor, £^ v , on the brane which is 
responsible for carrying on the brane the contribution of the bulk gravitational field. The cosmolog- 
ical evolution of such a brane universe has been investigated and effects such as a quadratic density 
term in the Friedmann equations have been found [6, 7, 8]. 

Lorentz symmetry is assumed to be the exact symmetry of nature [9]. However, many exotic 
theories such as canonical quantum gravity and string theory suggest that Lorentz invariance may 
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be broken at high energies [10, 11]. String theory predicts that we may live in a universe with 
non-commutative coordinate [12] leading to violation of Lorentz invariance as an explanation for the 
astrophysical anomalies such as the missing Greisen-Zatsepin-Kuzmin cutoff [13, 14, 15, 16]. A method 
of implementing the local Lorentz violation in a gravitational setting is to consider the existence of 
a tensor field with a non- vanishing expectation value and then couple it to gravity or matter fields. 
The simplest example of this method is to consider a single time-like vector field with fixed norm. 
A special case of this theory as a mechanism for Lorentz invariance violation has been introduced 
by Kostelecky and Samuel [17]. In a different context, a theory of gravity with a fixed norm vector 
in order to mimic the effects of dark matter has been proposed by Bekenstein [18]. Also, for the 
studies of vector fields in a cosmological setting without the fixed norm see [19, 20]. This vector field 
picks out a preferred frame at each point in space-time and any matter field coupled to this vector 
field will experience a violation of local Lorentz invariance [21]. However, in curved space-time there 
is no natural generalization of the notion of a constant vector field (V \ 1 N U = generically has no 
solutions); therefore we must allow that the vector field to have dynamic and fix the norm of it by 
choosing an appropriate action for the field. 

The effects of Lorentz violation have been also studied within the context of the brane-world 
scenarios [22]. In these models, the space-time violates four-dimensional Lorentz invariance globally 
and leads to apparent violations of Lorentz invariance from the brane observers point of view due 
to bulk gravity effects. These effects are restricted to the gravity sector of the effective theory while 
the measured Lorentz invariance of particle physics remains unaffected in these models [23]. In a 
similar way, Lorentz invariance violation has had to shed some light on the possibility of signals 
travelling along the extra dimension outside our universe [24]. In a different method a brane-world 
toy model in an inflating five-dimensional brane-world scenario has been introduced with violation 
of four-dimensional Lorentz invariance at an energy scale k [25]. In [26], the authors have studied 
Lorentz violation in a brane-world model by considering a vector field normal to the brane along the 
extra dimension. They showed that local Lorentz violation in the bulk allows for the construction 
of models in which the vacuum energy, the gravitational coupling and the cosmological term on the 
brane are variable. There also exist a relation between the maximal velocity in the bulk and the speed 
of light on the brane [27]. The main aim of this paper is to see how the shear parameter behaves 
in such brane-world scenario and whether this model isotropises at late time. To achieve this goal 
we study the local Lorentz invariance violation in the same setting and investigate the behavior of 
the observationally important parameters like shear, anisotropy, and the deceleration parameter in a 
Bianchi type I geometry. 



2 The setup 

Let us start by presenting the model used in our calculations. We only state the results and refer the 
reader to [26] for a detailed derivation of these results. 

As mentioned above, the brane-world model we consider here includes Lorentz violation in the 
bulk space along the extra dimension by generalizing the theory suggested by Jacobson and Mattingly 
[28, 29], where the quartic self-interaction term (N a VaN b )(N c X<'cNb) has been ignored [30]. We 
assume that N A is a vector field along the extra dimension which is making the associated frame a 
preferred one. We also consider that the theory which is consists of the vector field N A minimally 
coupled to gravity with the following action 

S = J d 5 X^& -L (^^ + C n) + Cm , (1) 

where k| = 8ttG^,, is the five-dimensional Ricci scalar, Cn is the vector field Lagrangian density 
and C m is the Lagrangian density for the matter fields. We take N A as a dynamical field for preserving 
the general covariance. The Lagrangian density for the vector field is given by 

C N = K AB CD V A N C V B N D + X(N A N A - e) , (2) 
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where A is a Lagrange multiplier and 

K AB CD = -Pig AB gcD - hs&l - faS&E, 



(3) 



here, f3i,i = 1,2,3, dimensionless parameters, e = —1 or e = +1 depending on whether the extra 
dimension is time-like or space-like, respectively. Using the metric of the bulk space as 



and 



dS 2 = g^(x a , y)dx^dx u + e0(x a , y)dy 2 , 



9 abN a N b = e, e 2 = 1, 



(4) 



(5) 



and assuming that the normal unit vector N A which is orthogonal to the hypersurfaces y = const. 
as, N = 4-, the gravitational field equations on the brane become as [26] 



G 



k 



3(e + a.i) 



3(e + a 3 ) 



9^uM - —— (K^KJ - KKfj, v ) - — — rg^K 



(3 + ai) 
3(e + a 4 ) K K a{S , 



2(3 + ai) 



&l(0l5 + \) + 3CK5 



(3 + ai) 



a 5 (5 + ai) <P, a (t) a 



2(3 + a 1 )^"</> : 



2ai 4a 5 _ 3(e + ai) 



where 



(3 + ai) (3 + ai) </> 2 (3 + ai) 

2e(fa + /3 2 + /3 3 



(6) 



ai = 2(/3i+/3 3 ), a 2 



3-2e(/3i+4/3 2 + /? 3 )' 



«3 



ai(3 + e-2/3 2 ) ai (6 + e + ai ) 

P2, «4 = - , a 5 = epi. 



6 " ' 6 

The first partial derivatives in terms of the extrinsic curvature of the brane is given by 



K, 



2<j> dy ' 



K A5 = 0, 



and 



c _ M 5 ) ajA N B 



(7) 



(8) 



(9) 



Using equation (6) and the Z 2 symmetry, we obtain the gravitational field equations on the brane as 



3(e + a±] 



(3 + ai) 



"1(05 + \) + 3tt5 

(3 + ai) 



a 5 (5 + ai) 0,a^° 



2a 



2(3 + ai)"^" 2 (3 + ai) <A ' (3 + ai) (f> 2 ' 

"e — 3(— ol\ + 4q 6 + 16«7 + 2a 4 ) 



1 <?W 



+ 



4a K 



where 



A 4 = ^A 5 + 4 



2 _ 4 

— Kr^ 



4(3 + ai)(l + eai) 2 
2e - 3(-2qi + 4a 6 



A 2 , 



4(3 + ai)(l + eai) 2 



A, 



(10) 

(11) 
(12) 



7T 



+ 



4(3 + ai)(l + eai) 2 
3(q6 + 807 + a 4 ) 
4(3 + ai)(l + eai) 2 



-(e + ai)r^ Q r" + (| + a 6 )rr^ - (| - a 7 )g flv T 2 + ^ g^T a pT al3 



gtiu^r, 



(13) 
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TP 

- 1 (IV 



\y=0, 



(14) 



and 

ai(l - 2a 2 ) - 2ea 2 , 1t .. 

«6 = g , (15) 

(e + ai)(a 2 + a 2 ) (a 4 - 3e - 4a 3 )(a 2 + 2a 2 ) (a 3 + 2a 4 ) 

" 7 = 3 + 9 18 ' (16) 

For latter convenience we chose ao = 2/3 2 , so that constants Qj (i = 2, 7) are defined as a function 

of «o and a\. Now, by considering N§ = 5$, e = +1 and ignoring the bulk matter, = 0, we 
obtain the gravitational field equations in four-dimensions as 



where 



Gfiv = -^49fMu + K4 ( V> + j^9^ur ) + K 5 vr Mi , - Sfj, v , (17) 

A4 = |a 5 + rA 2 , (18) 

2 4(3 - 4ao - on) 

K ' = 2(3 + a 1 )(3-4a -a 1 ) A ' (19) 

3 r (l-2a -ai) __ _ , (6 + Qi) „ a/ j 
4(3 + ai)(l + ai) (3 — 4ao — ai) p 12 
2(3 -ai) - (9 + ai)a ^ 2 
12(3-4a -ai) ^ ] ' (20) 



and 



£ _3(l + ai) f . 

We note that the model is different from the SMS model in two cases. The first one is the existence 
of the trace part of the brane energy-momentum tensor in the modified gravitational field equations 
on the brane. This trace part of the energy-momentum tensor vanishes when a,\ = 0. The second 
departure from the SMS model arises from definition of the fundamental quantities A4 and k\ which 
contain higher-dimensional modifications to the standard general relativity. 

The Codazzi equation also implies the conservation of the energy-momentum tensor of the matter 
on the brane 

V M r^ + a 2 V u T - (1 + 4a 2 )V,A = -2(e + ai )V v , (22) 

thus the brane energy- momentum tensor r fl u is not conserved. Moreover, the contracted Bianchi 
identities on the brane imply that the projected Weyl tensor should obey the constraint 

= -V„A 4 + k\ (V^ + gv„r) + k*VX„ + f V ^ + ^ ( V m^ " ^V„T" a ) . (23) 

For convenience we suppose that the bulk cosmological constant and brane tension is constant. Thus 
equations (22) and (23) become 

V M r^ = -a 2 V„T, (24) 

= -k\ (a 2 - V„r + «|V M 7r^. (25) 

Equations (17), (24) and equation (25) give a complete set of field equations for the brane gravitational 
field. In the next section, we study the cosmological consequence of anisotropic brane in the framework 
of such model. 
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3 Bianchi type I brane cosmology 



In the following we will investigate the influence of the Lorentz violating parameters on the anisotropic 
universe described by Bianchi type I geometry. Considering = 0, i = 1,2,3, we can reduce the 
model to the SMS brane-world model and compare our results to this model. 

The line-element of a Bianchi type I space-time, which is a generalization of the isotropic flat 
Friedmann-Robertson- Walker (FRW) metric, is described by 



ds 2 = -dt <3dt + ^2 ai{t) 2 dx l <g> dx\ 



(26) 



where ai(t),i = 1,2,3 are the expansion factors in different spatial directions. For later convenience 
we define the following variables 



i=i 



(27) 



Hi = 



i = 1,2,3, 



(28) 



(29) 



i=l 



In above equation, v is the volume scale factor, Hi,i = 1,2,3 are the directional Hubble parameters, 
and H is the mean Hubble parameter. The physical observable are the mean anisotropy parameter 
A, the deceleration parameter q and the shear scalar a 2 which are defined as 



/AH t 



3^V H 



(30) 



(31) 



VI 



l -[Y,H 2 -m 2 



(32) 



u=l 



where AHi = Hi — H . The sign of the deceleration parameter indicates how the universe expands. A 
positive sign for q corresponds to a decelerating universe whereas a negative sign indicates inflation. 
We note that A = for an isotropic expansion. 

We also assume that the confined matter source on the brane is a perfect fluid with equation of 
state p = (7 — l)p where 1 < 7 < 2. 

Using the variables (27)- (29), the Einstein gravitational field equation (17), the Bianchi identity 
(24) and the evolution equation for non-local dark radiation (25) take the form 



m + Y^H 2 = a 4 -| 
i=l z 



(3 7 -2) + ^(3 7 -4) 



K5(a 8 +3a 9 ) 2 3(1 + ai 



(3 + Ql) 



■U, (33) 



--(vHi) = A 4 -^ 
vdV l> 2 



(7 -2) + ^(3 7 -4) 



P + 



K^(a 8 - a 9 ) 2 (1 + ai] 



(3 + ai 



U, 



(34) 



5 



and 



p [l + a 2 (4- 3 7 )] + 3^Hp = 0, 



U + AHU 



«h(«2-^)(4-37)(3 + ai) 



«;|aio(3 + ai) 



(l + ai)[l + a 2 (4-37)] (1 + Qi)[l + a 2 (4 - 3 7 )] 

where 3-ff = | and constants as, ag and «io are defined as 

-4(1 + ai)(3 + ai) + 67(1 + ai)(4 + ai) - 37 2 (5 + ai)(a + «i) 



«8 



a 9 



16(3 + «i)(l + ai)(3 - 4a - aj) 
4(1 + ai)(3 + ai) - 6701(1 + ai) - 37 2 (a + «i)(3 - a{) 



16(3 + ai)(l + ai)(3 - 4a - ai) 
a i0 = 2 7 a 8 + (a 8 + a 9 )[l + a 2 (4 - 37)]. 
For 77^0 the solution of equation (35) is given by 



p = p Q V l+(4-3T)° 



«W"2-^)(4-3 7 )(3 + ai) T _ i ^ 



Also integrating equation (36) yields 
W = -3ctr 4 / 3 + 



^5Po a io(3 + ai) 



[4 + 4a 2 (4-3 7 ) -3 7 ](l + ai; 



- W l+(4-37)c«2 



(35) 
(36) 



(37) 



[4 + 4(4-3 7 )a 2 -6 7 ](l + ai 



. u l + (4-37)a 2 ^ 

(38) 

where c is a constant of integration. Using the non-local energy density and the evolution law of 
energy density, the gravitational field equations (33) and (34) become 



3H + J2 H i = A 4 - 



K4P0 



i=i 



2 

3(1 + a X/ 
(3 + ai) 



(3 7 -2) + ^(3 7 -4) 



, 1 + (4-37)q 2 



4p 2 o 



-27 



(a 8 + 3a 9 )v 1+ ( 4 - 3 ^ Q 2 



1,2,3, 



(39) 



f at 



2 

(1 + ai) 



(3 + ai) 

Summing equations (40) we find 
vdV ' 

+ 



(7 -2) + ^(3 7 -4) 
i = 1,2,3. 



w l + (4-3 7 ) Q2 + l5P0( a8 _ ag)?; l+(4-37)«2 



3A 4 



3«fpg 



3/?4 



(7 -2) + ^(3 7 -4) 



, l+(4-37)c 



-2 7 



3(1 + ai) , , 

(a 8 - ag)v T +^=^ r 2 - -A l ±u. 



2 v " u " (3 + ai) 

Now, substituting back equation (41) into equations (40) we obtain 



hi 

Hi = H -\ — -, 

v 



1,2,3, 



(40) 



(41) 



(42) 



with hi, i = 1,2,3 being constants of integration satisfying the consistency condition J2i=i hi = 0. The 
basic equation describing the dynamics of the anisotropic brane-world model with Lorentz violation 
in the bulk can be written as 



v = 3A 4 v + -^- 1/3 

(3 + ai) 



(l-7) + (4-37)a 2 (l-2 7 ) + (4-37)a 2 

+ /(a ,ai, 7 )t; i+(4-3 7 )«2 +T(a ,a 1 ,-y)v ^-^2 , (43) 
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where 



/(a ,O!i,7) = K 5 Ap 



, 18(1 + ai)[(2 - 7) + ^(4 - 3 7 )] - 7 (12q + 9<*i + a? + 4a ai) 



8(l + ai)(3 + ai)(3-4a -ai) 



J"(a ,ai,7) = «5Po 



3[4(3 + ai) + 37 2 (a + ai) - 67(2 + ai)] aio(3 - 4a - ai) 

16(3 + qi)(3-4q -Q!i) 4(1 + ai ) + 27(2a - a~i - 3) 



The general solution of equation (43) becomes 

I - , , 27c(l+«l) ,,„ . , (2-7)+2(4-3 7)a2 

t-to = /(3A 4 f 2 + — ^ r±V /3 + 2g(a ,a 1 ,j)v ^-^2 

J (6 + ai) 

2(l- 7 )+2(4-3 7 )a 2 . .„ 

+ 20(a o ,ai,7)u r+^^S + Cy^dv , (44) 

1 + (4 - 3 7 )a 2 
5(a °' ai ' 7) = (2-7) + 2(4-37)a 2 /(a °' ai ' 7) ' 

r , x 1 + (4 - 3 7 )a 2 

^ (a °' ai ' 7) = 2(l-7) + 2(4-37)a/ (a °' ai ' 7) ' 
where C is a constant of integration. The time variation of the physically important parameters 
described above in the exact parametric form, with v taken as a parameter, is given by 



where 



A = 3h 2 



3A 4 v 2 + 



27cil+aA 0/0 (2- 7 )+2(4-37)a 2 2(l- 7 )+2(4-3 7 ) a2 

; 1 i tly 2 ' 3 + 2g(a ,a 1 ,~ ( )v w-w* + 2Q(a , ai, j)v ^=^2 +C 



(3 + ai) 



-1 



,(45) 



3v 



(l_ 7 ) + (4-3 7 )a 2 



(l-2 7 ) + (4-3 7 )a 2 



3A 4 ^ + $+ai) v + /(«o,«i,7)^ 1+ ( 4 - 3 ^ Q 2 +7 , (a ,ai,7)u i+(4-3 7 )« 2 



(46) 



,-,,,„ x , (2- 7 )+2(4-3 7 )a 2 

3A 4 « 2 + ^ 2/3 + 2 5(«o, ai,7)u i+(4-37)o 2 + 2£(a , a u -y)v 



2(l- 7 )+2(4-3 7 )a 2 

l+(4-3 7 ) a2 _|_ Q 



r 27c(l+CKl) (2- 7 )+2(4-3 7 )ct 2 

<H = a 0i v l/3 exp{ / [3A 4 v 4 + — j- ^V /3 + 2 5 (a , 7)" 1+ < 4 " 3 ^ 



2(l- 7 )+2(4-3 7 )a 2 

+ 20(ao,ai,7)u i+(4-3 7)Q2 + C?; 2 ]~ 1/2 (iw}, i = 1,2,3, 



(47) 



a 2 = 



2v 2 ' 



(48) 



where h 2 = J2i=i hf- Also, the integration constants hi and C must satisfy the consistency condition 
h 2 = %C. As one can see from equations (45)-(47) the behavior of these physical parameters depends 
on the Lorentz violating parameters and the equation of state of the cosmological fluid. 
For a well-defined theory, the following constraints put on the parameters ft as [30] 

ft >o, 



(ft + ft + ft)/ft <i, 
(ft + ft + ftVA > 0, 
ft + ft < 0. 
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constraints for G > G constraints of ref [30] integer values of «o and a.\ 

which satisfy both of two constraints 



a > n 


3 


— 4ao — «i > 


nnci "H VP 
JJVJol LI V 


rvn > D 
«0 + Qi > 


(Tn — D (T1 — D 1 2 






3 + ai < 


positive 


a > 


a = {4,5,6, ...} 




3 


— 4qo — ol\ < 




«o + a i > 


ol\ = {-a , -a + 1, -a + 2..., -4} 


A < 


3 


3 + ai < 
— 4qo — ai > 


positive 


«o > 
+ a i > 


no values 






3 + ai > 


positive 


a > 


ao = 0, a.\ > 3 




3 


— 4ao — «i < 




«0 + CKi > 


ao = 1, ot\ > —1 



a = {2,3,4,...}, ai > -3 



Table 1: The possible values of Lorentz violating parameters in a brane- world universe. 



The first condition results from the need for a positive-definite Hamiltonian for the perturbations; 
the next two from demanding the subluminal and non-tachyonic propagation of the spin-0 field; and 
the last condition from insisting subluminal propagation of the spin-2 field. Together these conditions 
imply that 

Pi + h + fa > 0, 



which in terms of ao and ai we have 



ft >o, 

ao + a i > 0, 
a > 0. 



The above conditions together with constraints require for a positive effective Newtonian constant 
(19), lead to the particular values of ao and a±. We have summarized these results in table 1. 

The singular state at t = is characterized by the condition v(0) = 0. The value of the anisotropy 
parameter for t = depends on the Lorentz violating parameters and the equation of the state. Hence 
for 7 = 1 equation (45) reduces to 

3A 4 ^ 2 + ( 3 ( ++ ) l} ^ + 2g(a , a u l)v ^ + 2£(a , «i, l)v ^ +C 

the singular behavior of the dust filled brane universe depends on 02, so that from above equation it 
follows 

lim A(v) = 0, -1< a 2 < 0, (50) 

v— >0 

and 

36/i 2 

S A( » ) = fcfpS + 12C7 ' O2 = ' (51) 

and otherwise 

3/) 2 

UmA(v) = ^. (52) 
v^O C 

In the absence of the Lorentz violating parameters the dust filled brane universes are born in an 
anisotropic state [32] , while in the brane universe with the Lorentz invariance violation the early time 
behavior of the brane universe is sensitive to the values of a 2 , so that it is possible to admitting of 
both the isotropic and anisotropic state. Therefore the singular state of the dust filled Bianchi type I 
brane universe in our model is different from the case of brane cosmological models without the bulk 
Lorentz violation. 
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For the case of the radiation anisotropic brane universe, p = ^p and 7=3, equation (45) reduces 

to 

3/? 2 

A = ™ (53) 

3A 4 t; 2 + ^IT^ vl + 2g(a , a u 4/3)«l + 2G(a , a u 4/3)v~l + C 

in this case the singular behavior of the anisotropy parameter is independent of «2, so that from 
equation (53), we have 

lim A(v) = 0. (54) 

V— ¥0 

In figure 1, we have plotted the anisotropy parameter for cto = and a\ = 2. Since 02 expressed 
in terms of Qo and a.\ becomes 02 = 3 _" 4 Q ^" 1 Qi , therefore «2 = 2 and from equation (52) for 7 = 1 

we have A(0) = and the anisotropic brane universe starts its evolution from an anisotropic state 
to an isotropic de Sitter inflationary phase at late time. For 7 = 4/3, the late time behavior of 
anisotropy parameter is isotropic with A(0) = 0. The behavior of the deceleration parameter of the 
Bianchi type I geometry is illustrated, for different values of 7, «o = and a\ = 2, in figure 2. In the 
initial stage the behavior of the Bianchi type I brane universe is non-inflationary, but brane universe 
ends in an accelerated expanding stage at late time. The time variation of the shear parameter as a 
function of time is represented in figure 3. The shear is a decreasing function of time and in the limit 
of large time corresponding to the isotropic limit, o 2 — > 0. In the limit of small time the shear has a 
singular behavior tending to infinity, a 2 — > 00. 

For fa = and fi\ = — /?3, equations (45)-(47) give the general solutions of the physical parameters 
on the volume scale factor for the generalized RS II model [33] 



aoiV 1 ^ 3 exp 



hi J (?,k A v 4 + ?,K 2 p v 4 -^ + ^4 PoV 4 ~ 21 + Cv 2 ^j 1 dv 



1,2,3, (55) 



_ 36A 4 ^ 2 + 18(2 - 7 )k 2 pqv 2 -^ + 3(1 - 7)^Po^ 2 " 27 ^ 

Q 12A 4 V 2 + 12/^^2-7 + 4^2-27 + 4C ' ^ 

A = 3/i 2 (3A 4 v 2 + 3^p v 2 " 7 + ^4p 2 v 2 - 2 ~< + C^j 1 , (57) 



where h 2 = X)f=i ^1 an d C is a constant of integration. The behavior of the deceleration parameter 
of the Bianchi type I geometry is illustrated, for different values of 7, in figure 4. In the initial stage 
the behavior of the Bianchi type I brane universe is non-inflationary, but brane universe ends in an 
accelerated expanding stage at late time. 

In figure 5, we have also plotted the anisotropy parameter for different values of 7. The behavior 
of the anisotropy parameter shows that at high densities the brane universe starts its evolutions from 
an isotropic state with A(Q) = for 7 = 4/3, and ends in an isotropic de Sitter inflationary phase at 
late time. An important difference between the anisotropic homogeneous brane-world cosmological 
models and the standard general relativity is that brane universes are born in an isotropic state. 

The time variation of the shear parameter as a function of time is represented in figure 6. The 
shear is a decreasing function of time and in the limit of large time corresponding to the isotropic 
limit, a 2 — > 0. In the limit of small time the shear has a singular behavior tending to infinity, a 2 — > 00. 
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A(t) 




Figure 1: The anisotropy parameter of the Bianchi type I brane geometry, with the bulk Lorentz violation ao = 0, 
ai = 2, for 7 = 4/3 (solid curve) and for 7 = 1 (dashed curve). 



q(t) 




Figure 2: The deceleration parameter of the Bianchi type I brane geometry, with the bulk Lorentz violation ao = 0, 
ai = 2, for 7 = 4/3 (solid curve) and for 7 = 1 (dashed curve). 

4 Conclusions 

In this letter, we have dealt with the dark energy problem in the context of the brane- world scenario 
with the bulk Lorentz violation, introduced by specifying a preferred frame through the introduction 
of a dynamical vector field normal to the brane. In our model, due to the local Lorentz violation 
in the bulk, the Friedmann equations on the anisotropic brane have been modified by the Lorentz 
violating parameters /3j,i = 1,2,3 and the equation of state of the matter. Therefore, for a fixed 
value of 7, the behavior of an anisotropic brane universe is controlled by the Lorentz violating pa- 
rameters. The behavior of the anisotropy parameter shows that the brane universe evolves from an 
isotropic/anisotropic state to an isotropic state which has entered an accelerated expanding phase. 
The early time behavior of anisotropic parameter is different from those obtained in RS type II 
brane- world models without the bulk Lorentz violation [33]. 
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Figure 3: The shear parameter of the Bianchi type I brane universe, with the bulk Lorentz violation ao = 0, ai = 2, 
for 7 = 4/3 (solid curve) and for 7 = 1 (dashed curve). 
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Figure 4: The deceleration parameter of the Bianchi type I brane universe with confined perfect cosmological fluid for 
7 = 4/3 (solid curve) and 7 = 1 (dashed curve). 
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Figure 5: The anisotropy parameter of the Bianchi type I brane universe with confined perfect cosmological fluid for 
7 = 4/3 (solid curve) and 7 = 1 (dashed curve). 
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Figure 6: The shear parameter of the Bianchi type I brane universe with confined perfect cosmological fluid for 7 = 4/3 
(solid curve) and 7 = 1 (dashed curve). 
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